Two-dimensional Potts antiferromagnets with a phase transition at arbitrarily large q 
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We exhibit infinite families of two-dimensional lattices (some of which are triangulations or quad- 
rangulations of the plane) on which the q-state Potts antiferromagnet has a finite-temperature phase 
transition at arbitrarily large values of q. This result is proven rigorously using a Peierls argument. 
Additional numerical data are obtained using transfer matrices, Monte Carlo simulation, and a 
high-precision graph-theoretic method. 

PACS numbers: 05.50.+q, ll.10.Kk, 64.60.Cn, 64.60.De 



X 



The q-state Potts model [1, 2] plays an important 
role in the theory of critical phenomena, especially in 
two dimensions [3-5], and has applications to various 
condensed- matter systems [2] . Ferromagnetic Potts mod- 
els are by now fairly well understood, thanks to universal- 
ity; but the behavior of antiferromagnetic Potts models 
depends strongly on the microscopic lattice structure, so 
that many basic questions about the phase diagram and 
critical exponents must be investigated case-by-case. 

One expects that for each lattice C there exists a value 
q c {£) [possibly noninteger] such that for q > q c (C) the 
model has exponential decay of correlations at all tem- 
peratures including zero, while for q = q c (C) the model 
has a zero-temperature critical point. The first task, for 
any lattice, is thus to determine q c . 

Some two-dimensional (2D) antiferromagnetic mod- 
els at zero temperature can be mapped exactly onto a 
"height" model (in general vector- valued) [6, 7]. Since 
the height model must either be in a "smooth" (or- 
dered) or "rough" (massless) phase, the corresponding 
zero-temperature spin model must either be ordered or 
critical, never disordered. Experience tells us that the 
most common case is criticality [8]. The long-distance 
behavior is then that of a massless Gaussian with some 
(a priori unknown) "stiffness matrix" K > 0. The crit- 
ical operators can be identified via the height mapping, 
and the corresponding critical exponents can be pre- 
dicted in terms of K. Height representations thus provide 
a means for recovering a sort of universality for some (but 
not all) antiferromagnetic models and for understanding 
their critical behavior in terms of conformal field theory. 

In particular, when the y-state zero-temperature Potts 
antiferromagnet (AF) on a 2D lattice £ admits a height 
representation, one ordinarily expects that q = q c (£)- 



This prediction is confirmed in most heretofore-studied 
cases: 3-state square-lattice [6, 9, 12, 13], 3-state kagome 
[14, 15], 4-state triangular [16], and 4-state on the line 
graph of the square lattice [15, 17]. Until recently the 
only known exception was the triangular Ising AF [18]. 

Kotecky, Salas and Sokal (KSS) [10] have observed 
that the height mapping employed for the 3-state Potts 
AF on the square lattice [6] carries over unchanged to any 
plane quadrangulation; and Moore and Newman [16] ob- 
served that the height mapping employed for the 4-state 
Potts AF on the triangular lattice carries over unchanged 
to any Eulerian plane triangulation (a graph is called 
Eulerian if all vertices have even degree). One therefore 
expects naively that q c = 3 for every (periodic) plane 
quadrangulation, and that q c = 4 for every (periodic) 
Eulerian plane triangulation. 

Surprisingly, these predictions are falsel KSS [10] 
proved rigorously that the 3-state AF on the diced lat- 
tice (which is a quadrangulation) has a phase transi- 
tion at finite temperature (see also [20]); numerical esti- 
mates from transfer matrices yield g c (diced) ss 3.45 [21]. 
Likewise, we recently [11] provided analytic arguments 
(falling short, however, of a rigorous proof) that on any 
Eulerian plane triangulation in which one sublattice con- 
sists entirely of vertices of degree 4, the 4-state AF has a 
phase transition at finite temperature, so that q c > 4. We 
also presented transfer-matrix and Monte Carlo data con- 
firming these predictions for the union-jack and bisected 
hexagonal lattices, leading to the estimates g c (UJ) sa 4.33 
and 9 C (BH) 5.40. 

These results suggest the obvious question: How large 
can q c be on a plane quadrangulation (resp. Eulerian 
plane triangulation)? The answers are clearly larger than 
3 or 4, respectively — but how much larger? 



Figure 1. Unit cells of the lattices G n and H„ for n — 3. 

In this Letter we shall give a rigorous proof of the unex- 
pected answer: we exhibit infinite classes of plane quad- 
rangulations and Eulcrian plane triangulations on which 
q c can take arbitrarily large values. We shall also comple- 
ment this rigorous proof with detailed quantitative data 
from transfer matrices, Monte Carlo simulations, and a 
powerful graph-theoretic approach developed recently by 
Jacobsen and Scullard [22]. 

The lattices G n and H n . Let G n be obtained from 
from the square (SQ) lattice by replacing each edge with 
n two-edge paths in parallel; and let H n be obtained from 
G n by connecting each group of n "new" vertices with an 
(n— l)-edge path (see Fig. 1). Resumming over the spins 
on the "new" vertices [23], it is easy to show that the 
q-state Potts model on G„ or H n with nearest-neighbor 
coupling v — e J — 1 is equivalent to a SQ-lattice Potts 
model with a suitable coupling V e g(q,v) [24]; moreover, 
for q > 2 (rcsp. q > 3) an AF model (—1 < v < 0) on G n 
(resp. H n ) maps onto a ferromagnetic model (v e g > 0) 
on the SQ lattice. Concretely, for the zero-temperature 
AF (v = — 1) we have 

»£-ta,-l) = (^i)"-l (1) 

Setting v e s equal to the SQ-lattice ferromagnetic critical 
point v c (SQ) = y/q [3, 25], we obtain q c for G n and H n ; 
they have the large-n asymptotic behavior 

In 

q c (G n ) « q c (H n ) » + 0((n/\ogn)^) (3) 

where W(x) w \ogx — log log x + o(l) is the Lambert 
W function [26]. We have thus exhibited two infinite 
families of periodic planar lattices on which the Potts 
AF has arbitrarily large q c as n — > oo [27]. These lattices 
are not triangulations or quadrangulations, but they can 
be modified to be such and retain the phase transition, 
as we now show. 

The modified lattices. Starting from G n or H n , insert 
a new vertex into each octagonal face and connect it ci- 
ther to the four surrounding vertices of the original SQ 
lattice, to the four "new" vertices, or to all eight vertices; 
call these modifications ', " and "', respectively In par- 



Figure 2. Unit cells of the lattices G n , G'n and Hn for n = 3. 

ticular, G' n and G" are quadrangulations, and H^' is an 
Eulerian triangulation (Fig. 2). 

If we integrate out the spins at the vertices placed into 
the octagonal faces, we obtain the model on G n or H n 
perturbed by a 4-spin or 8-spin interaction. When q is 
large, this interaction is weak (of order 1/q) because its 
Boltzmann weight is bounded between a maximum value 
of q and a minimum value of q — 4 or q — 8. We therefore 
expect that the new edges will have a negligible effect 
on the phase transition when q is large, and that all the 
modified lattices will have q c (n) whose large- n behavior 
is essentially identical to Eq. (3). Let us now sketch a 
rigorous proof [28] of this assertion. 

Proof of phase transition. Recall first how one proves, 
using the Peierls argument, the existence of ferromag- 
netic long-range order (FLRO) at low temperature in the 
q-state Potts ferromagnet on the SQ lattice. The Peierls 
contours are defined as the connected components of the 
union of all bonds on the dual SQ lattice that separate 
unequal spins. A Peierls contour 7 of length I7I and cyclo- 
matic number 0(7) comes with a weight that is bounded 
above by (<jr— l) c ^ 7 ^(l + w)"' 7 ' : here (q— 1) c (t) is a bound 
on the number of colorings of the SQ lattice consistent 
with the contour 7. Further, on the SQ lattice we have 
c {l) < ItI/2, and the number of contours of length n 
surrounding a fixed site can be bounded by (n/2)16". 
Standard Peierls reasoning then shows that for any pair 
of sites x, y one has 

00 

Prob(<7 K ^ o v ) < J2(n/2)16 n (q-l) n / 2 (l+v)- n , (4) 

which is < 5/16 whenever 1 + v > 32^/q. This proves 
FLRO (the constant 32 is of course suboptimal). The 
foregoing argument is valid for fixed boundary conditions 
(e.g., a = 1) in the plane, but with suitable modifications 
it can also be carried out for periodic boundary condi- 
tions (i.e., on a torus). 

Let us now consider the Potts antiferromagnet on one 
of the six modified lattices G' n , . . . , H"' . Since our goal 
is to show FLRO on the SQ sublattice, we define Peierls 
contours exactly as we did for the SQ-lattice ferromag- 
net, ignoring the spin values at all other sites. Although 
we no longer have any simple explicit formula for the 
contour weights, it is nevertheless possible to prove an 
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Figure 3. The union-jack (UJ) lattice. 



upper bound on the probability of occurrence of a con- 
tour 7 by using the technique of reflection positivity 
and chessboard estimates [29]. Without going into de- 
tails of the needed adaptations of this standard tech- 
nique for our case (see [28]), we mention only that the 
final bound on the probability of occurrence of a con- 
tour 7 is (n,y/q — l)' 7 ' , where k is the probability that 
the spins on the SQ sublattice follow a fixed checker- 
board pattern (say, 1 on the even sublattice and 2 on the 
odd sublattice) raised to the power 1/volume. This lat- 
ter probability is easy to bound explicitly, yielding k < 
[l+v c g(q, v)]^ 1 [q/ (q— 5)] 1 / 2 , where v ff (q, v) is the one for 
the corresponding unmodified lattice G„ or H„ . This im- 
plies that, for all the lattices G' n , . . . , if"', there is FLRO 
on the SQ sublattice whenever 6 < q < q c (G n ) — 0(1) 
[cf. Eq. (3)] and v is close to —1 (low temperature). 

Let us also remark that the lattice G'{ is covered by 
the general theory of [20], where it is proven that q c > 3; 
moreover, a minor modification proves the same result 
for G£ for all n > 2. 

Data for lattices G' n and H' n . The lattices G' n and 
H' n for all n can be reduced to the union-jack (UJ) lat- 
tice (Fig. 3) with v x = v and up = a suitable v c g(q 1 v) [cf. 
Eqns. (l)/(2) when v = —1]; of course the same reduc- 
tion holds for G n and H n by setting v x = 0. We obtained 
high-precision estimates of the phase boundary of the 
UJ model in the (w x , i>rj)-plane by using the Jacobsen- 
Scullard (JS) method [22] with untwisted square bases of 
size up to 6 x 6 (216 edges) [30]. We checked these re- 
sults for q = 4, 8, 16, 32 by Monte Carlo simulations using 
a cluster algorithm [31]. The estimated phase boundaries 
from both methods are shown in Fig. 4a, along with the 
numerical estimates of v\j c at v x = —1. The result- 
ing estimates of n c (q) [the inverse of the function q c (n)] 
from Eqns. (l)/(2) are shown in Fig. 4b, where they are 
compared with the predicted large-g asymptote n c (q) 
iqlogq from Eq. (3). The functions q c (n) divided by 
their large-n asymptote 2n/W(2n) are plotted in Fig. 4c. 
Note that q c (G' n ) > q c (G n ) and q c {H' n ) > q c (H n ), in ac- 
cordance with the intuitive idea that the AF edges asso- 
ciated to the modification ' enhance the ferromagnetic 
ordering on the SQ sublattice [32]. 
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Figure 4. (a) Estimated phase boundaries for q = 4, 8, 16, 32 
from the Jacobsen-Scullard method (blue curve and numer- 
ical values of vn iC at v x ~ —1) and Monte Carlo simula- 
tions (black points), (b) Estimates of n c (q), and their large-g 
asymptote from Eq. (3). (c) Plots of q c (n) divided by their 
large-n asymptote 2n/W(2n). 



Data for lattices G" and H"' . We studied the lattices 
G" and H% for n = 1, 2, 4, 8, 16, 32, 64 (note that G'{ = 
SQ [6] and H'{' = UJ [11]) at v = — 1, using transfer 
matrices with cylindrical boundary conditions on widths 
L = 1, 2, 3, 4 unit cells (Fig. 5). The computational com- 
plexity is linear in n. We estimated the location of the 
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Figure 5. The lattice (rotated 45° from Fig. 2) with 
L — 2; a unit cell is shown in red, and the transfer direction 
is indicated with an arrow. 

phase transition (which we expect to be first-order when- 
ever q c > 4) using the crossings of the energies El (q) [33] : 
the results are shown in Fig. 6a. 

For H%' we checked these results by Monte Carlo: for 
three integer values of q below the estimated q c we sim- 
ulated the model at finite temperature and estimated 
the transition point v c (q); we then performed linear and 
quadratic extrapolations to locate the point q c where 
v c = — 1. The results are shown in Fig. 6a, b and agree 
well with the transfer-matrix estimates. For q > 8 the 
specific heat diverges at the transition point like L~ 2 , 
in agreement with the finite-size-scaling prediction for a 
first-order transition; for 4 < q < 8 the transition is 
presumably also first-order but with a large correlation 
length £, so that we are unable to observe the true L ^> £ 
asymptotic behavior. 
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